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ABSTRACT 

We study the problem of detecting, and infering astrophysical information from, grav- 
itational waves from a pulsating neutron star. We show that the fluid / and p-modes, 
as well as the gravitational-wave u>-modes may be detectable from sources in our own 
galaxy, and investigate how accurately the frequencies and damping rates of these 
modes can be infered from a noisy gravitational-wave data stream. Based on the 
conclusions of this discussion we propose a strategy for revealing the supranuclear 
equation of state using the neutron star fingerprints: the observed frequencies of an 
/ and a p-mode. We also discuss how well the source can be located in the sky using 
observations with several detectors. 
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1 INTRODUCTION 



useful constraints on the equation of state ( Andersson & 



As we appiuach Ihu next millennium Llieie is a fucussed 

worldwide effort to construct devices that will enable the 
first undisputed detection of gravitational waves. A network 
of large-scale ground-based laser-interferometer detectors 
(L1GO, VIRGO, GEO600, TAMA300) is under construc- 
tion, while the sensitivity of the several resonant mass detec- 
tors that are already in operation continues to be improved. 
At the present time it seems likely that the gravitational- 
wave window to the Universe will be opened within the next 
five to ten years, and that gravitational- wave astronomy will 
finally become a reality. 

An integral part in this effort is played by theoretical 
modelling of the expected sources. Theorists are presently 
rackingi thoir braine to think of varioue eourcoe of gravita — 

tional that may ho nhgnrirahlp nnro tlio new ultraacn- 

sitive detectors operate at their optimum level, and of any 
piece of information one may be able to extract from such 
observations. One of the most challenging goals that can (at 
least, in principle) be achieved via gravitational wave detec — 
tion is the d£temunatkta of the equation " f state of mattpr 



at supranuclear densities. 

We have recently argued that observed gravitational 
waves from the various nonradial pulsation modes of a neu- 
tron star can be used to infer both the mass and the radius 
of the star with surprisingly good accuracy, and thus put 



Kokkotas 1996; Andersson & Kokkotas 1998). This conclu 



sion was, however, based on an "ideal" detection situation 
that ignored the various uncertainties associated with the 
analysis of a noisy data-stream. In the present paper we an- 
alyze thoroughly this idea by incorporating all possible sta- 
tistical errors that might arise in estimating the parameters 
associated with the neutron star equation of state. 

The spectrum of a pulsating relativistic star is known to 
be tremendously rich, but most of the associated pulsation 
modes are of little relevance for gravitational-wave detec- 
tion. From the gravitational-wave point of view one would 
expect the most important modes to be the fundamental (/) 
mode of fluid oscillation, the first few pre ssure (p) modes and 
the f irst gravitational- wave (w) modes (Kokkotas & Schutz 



1992 ). For details on the theory of relativistic stellar pulsa- 
t ion we refer the reader to a recent review article by one of us 
(Kokkotas 1997). That the bulk of the energy from an oscil- 
lating neutron star is, indeed, radiated through the se modes 
has been demonstrated by numerical experiments (Allen, et 
al. 1998j; [Allen, et al. 1999|). 



Recently, two of us (Andersson & Kokkotas 1998) pro- 



vided data for the relevant pulsation modes of several stel- 
lar models for each of twelve proposed realistic equations 
of state, thus extending earlier results of Lindblom and De- 
tweiler ( 1983 ) . This data was then used to obtain empirical 
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relations between the observables (frequency and damping 
time) o f the /- and the m-modes and the stellar parameters 



long enough (many dynamical timescales) that several os- 
cillation m odes could potentially be identified (Baumgarte 



(mass fnd radius). It was shown that these relations could et al 1996). Also, stellar oscillations can be excited by the 



be used to infer both the radius and the mass of the star 
(typically with an error smaller than 10%), i.e., to take the 
fingerprints of the star. It was also pointed out that, since 
no general empirical relations could be inferred for the p- 
modes, they could prove important for deducing the actual 
equation of state once the radius and the mass of the star is 
known. The empirical relations, including the relevant sta- 
tistical errors, are listed in Appendix 1. 

The proposed strategy can potentially be of great im- 
portance to gravitational-wave astronomy, since most stars 
are expected to oscillate nonradially. The evidence for this is 
compelling: Many p-modes (as well as possible (/-modes and 
r-modes) have been observed in the sun and there are strong 
indications that similar modes are excited also in more dis- 
tant stars. In principle, one would expect the modes of a 
star to be excited in any dynamical scenario that leads to 
significant asymmetries. Still, one can only hope to observe 
gravitational waves from the most compact stars. Hence, our 
attent ion is restricted t o neutron stars (or possible strange 
stars ( Alcock et al 1986| ), if they exist). Furthermore, to lead 
to detectable gravitational waves the modes must be excited 
to quite large amplitudes, which means that only the most 
violent processes are of interest. 

There are several scenarios in which the various pul- 
sation modes may be excited to an interesting level: (1) A 
supernova explosion is expected to form a wildly pulsating 
neutron star that emits gravitational waves. The current es- 
timates for the energy radiated as gravitational waves from 
supernovae is rather pessimistic, suggesting a total release of 
the equivalent to 10 -6 M c 2 , or so. However, this may be a be unstable (Andersson 1998; Friedman & Morsink 199£ 
serious underestimate if the gravitational collapse in which 
the neutron star is formed is strongly non-spherical. Opti- 
mistic estimates suggest that as much as 10 -2 Mqc? may 
be released in extreme events. (2) Another potential exci- 
tation mechanism for stellar pulsation is a starquake, e.g., 
associated with a pulsar glitch. The typical energy released 
in this process may be of the order of the maximum me- 
chanical energy that can be stored in the crust, estimated 



tidal fields of the two stars during the inspiral p hase that 
preceeds the merger ( [Kokkotas fe Schafer 1995] ). (4) The 
star may undergo a dramatic phase-transition that leads to 
a mini-collapse. This would be the result of a sudden soften- 
ing of the equation of state (for example, associated with the 
formation of a condensate consisting of pions or kaons). A 
phase-transition could lead to a sudden contraction during 
which a considerable part of the stars gravitational bind- 
ing energy would be released, and it seems inevitable that 
part of this energy would be channeled into pulsations of the 
remnant. Large amounts of energy that could be released in 
the most extreme of these scenarios: A contraction of (say) 
10% can easily lead to the release of 1O _2 M0C 2 . Transforma- 
tion of a neutron star into a strange star is likely to induce 
pulsations in a similar fashion. 

It is reasonable to assume that the bulk of the total 
energy of the oscillation is released through a few of the stars 
quadrupole pulsation modes in all these scenarios. We will 
assume that this is the case and assess the likelihood that 
the associated gravitational waves will be detected. Having 
done this we discuss the inverse problem, and investigate 
how accurately the neutron star parameters can be inferred 
from the gravitational wave data. 

Before we proceed with our main analysis it is worth- 
while making one further comment. In this study we ne- 
glect the effects of rotation on the pulsation modes. This 
is not because rotation plays an insignificant role. On the 
contrary, we expect rotation to be highly relevant in many 
cases. The re cent discovery th a t the so-called r-modes may 



and lead to rapid spin down of a young neutron star that i s 
born rapidly rotating (Andersson, Kokkotas & Schutz 1999), 



to be a t the level of 10 ~ 9 - 10~ 7 A/^c 2 (Blaes 1997 ; Mock & 



highlights the fact that it is not sufficient to consider only 
non-rotating stars. Still, as far as most pulsation modes are 
concerned, one would expect rotation to have a significant 
effect only for neutron stars with very short period, and the 
present study may well be reasonable for stars with periods 
longer than, say, 20 ms. Furthermore, it has been argued 
that neutron sta rs are typically born slowly rotating ( fgpruit 



Joss 19^7 ). This is also an interesting possibility consider- fc Phinney 1998 ). If that is the typical case, then our results 



ing the recent conclusion that the soft-gamma repeaters are 
likely to be so-c alled magnetars, neutron star s with extreme 
magnetic fields (Duncan & Thompson 1992), that undergo 
frequent starquakes. It seems very likely that some pulsation 
modes are excited by the rather dramatic events that lead to 
the most energ etic bursts seen from these systems. Indeed, 
Duncan (1992) has recently argued that toriodal modes in 
the crust should be excited. If modes are excited in these 
systems, an indication of the energy released in the most 
powerful bursts is the 1O _9 M0C 2 estimated for the March 
5 1979 burst in SGR 0526-66. The maximum energy should 
certainly not exceed the t otal supply in the magneti c field 
~ 10" 6 (B/10 15 G) 2 M Q c 2 (jDuncan & Thompson 19921). The 



possibility that a burst from a soft gamma-ray repeater may 
have a gravitational- wave analogue is very exciting. (3) The 
coalescence of two neutron stars at the end of binary inspiral 
may form a pulsating remnant. It is, of course, most likely, 
that a black hole is formed when two neutron stars coalesce, 
but even in that case the eventual collapse may be halted 



could be relevant for most newly born neutron stars. A more 
pragmatic reason for not including rotation in the present 
study is that detailed data for modes of rotating neutron 
stars is not yet available. Once such data has been com- 
puted the present study should be extended to incorporate 
rotational effects. 

Having pointed out this caveat, we are prepared to pro- 
ceed with the discussion of the present results. The rest of 
the paper is organized as fo llows . In Section 2 we basically re- 
peat the analysis of Finn (1992), in a slightly different form 
in order to compute the measurement errors of frequency 
and damping time of a gravitational wave that is emitted 
from a pulsating neutron star, and extend this analysis, so 
as to compute the accuracy by which one can estimate var- 
ious parameters of the star. In Section 3 we transform our 
results into a form that shows clearly how plausible it is 
to determine the equation of state of a neutron star by an- 
alyzing the gravitational wave data. In Section 4 we dis- 
cuss how detection of mode signals can be used to reveal 
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the nuclear equation of state. Section 5 is a brief discus- 
sion of how accurately one can expect to be able to locate 
the source in the sky. The final section presents our con- 
clusions. Appendix 1 contains empirical relations for oscil- 
lation frequency and damping rate of the / and ui-modes 
in terms of the stellar parameters (mass an d radius), de- 
duced from twelve realistic equations of state (Andersson & 



Kokkotis 1998). Appendix 2 lists the elements of the Fisher 



and covariance matrices used in the statistical analysis of 
the mode-signals. 



2 STATISTICAL ANALYSIS OF OBSERVED 
MODE-SIGNALS 

Suppose that one tries to detect the gravitational waves as- 
sociated with the stellar pulsation modes that are excited 
when (say) a neutron star forms after a supernova explo- 
sion. Since all modes are relatively short lived, the detection 
situat ion is similar t o that for a perturbed rotating black 
hole (Echeverrial989; Finnl992). For each individual mode 



(1) 



the signal is expected to have the following form: 

|0 fort<T, 
( ' ~ \ Ae- (t - T)/T sin[27r/(t - T)] for t > T. 

Here, A is the initial amplitude of the signal, T is its arrival 
time, and / and r are the frequency and damping time of 
the oscillation, respectively. Since the violent formation of a 
neutron star is a very complicated event, the above form of 
the waves becomes realistic only at the late stages when the 
remnant is settling down and its pulsations can be accurately 
described as a superposition of the various modes, either 
fluid or spacetime ones, that have been excited. At earlier 
times (t < T) the waves are expected to have a random 
character that is completely uncorrelated with the intrinsic 
noise of an earth-bound detector. This partly justifies our 
simplification of setting the waveform equal to zero for t < 
T. 

The energy flux F carried by any weak gravitational 
wave h is given by 



-M 



i, • ( 2 ) 

where c is the speed of light and G is Newton's gravitational 
constant. Thus, when gravitational waves emitted from a 
pulsating neutron star hit such a detector on Earth, their 



initial amplitude will be (Thorne 1987; Schutz 1997) 



A~ 2.4 



10" 



Etr 



io- 



/ lOkpc 



6 M c ; 
1kHz 



1/2 



)( 



lms 



(3) 



where _E gw is the energy released through the mode and 
r is the distance between detector and source. Again, / is 
the frequency of gravitational waves and r is the damping 
time: the rate at which the amplitude of the mode decays 
as energy is carried away from the star. 

In order to reveal this kind of signal from the noisy 
output of a detector one could use templates of the same 
form as the expected signal (so calle d mat ched filtering). 
Following the analysis of Echeverria ( 198£ ) the signal-to- 
noise ratio is found to be 



(I) 



N = P 



A 2 t 

1 + 4Q 2 2S n 



with 



Q = t/t , 



(4) 



(5) 



being the quality factor of the oscillation, and S n the spec- 
tral density of the detector (assumed to be constant over the 
bandwidth of the signal) . 

In Eq. (JiJ) , we used the following definition for the scalar 
product between two functions: 



{hi | h») = 



1 



Sn(\ f I) \ 

dthi(i)h a (t) 



(6) 



To compute the accuracy by which the parameters of the 
signal can be determined, we first define the dimensionless 
parameters e, -q, £, £ as 



foe 

To 77 



f - fo 
T — To 

T-T 



(7) 



where / , r , T , Ao are the true values of the four quantities 
/, r, T, A. These new parameters are simply relative devia- 
tions of the measured quantities from their true values. One 
can then construct the Fisher information matrix Tij and by 
inverting it, obtain all possible information about the mea- 
surement accuracy of each parameter of the signal, and the 
correlations between the errors of the parameters. The com- 
ponents of the symmetric Fisher matrix, which is defined 
by 



dh , dh 



(8) 



where 9i — (e, r/, are the parameters of the signal, are 
listed in Appendix 2. 

The inverse of the Fisher matrix, Ej 3 - = F" 1 , the so- 
called covariance matrix, is the most important quantity 
from the experimental point of view. Its components, that 
are directly related with the measurement errors of the pa- 
rameters, are also listed in Appendix 2. The purpose of the 
present analysis is to try to identify the nuclear equation 
of state from supposedly detected mode data. The relevant 
parameters for this analysis are the frequency / and damp- 
ing time r. Hence, we only need E e< ; and ~S vrl — the squares 
of the relative errors of / and r respectively — from Ap- 
pendix 2. We will also discuss the possibility that the signal 
will be seen by several detectors. If that is the case, we can 
use the time of arrival T to locate the position of the source 
in the sky. For this discussion we need E^, see Section 5. 

The estimated errors of the actual measurements should 
be taken into account along with the statistical errors in our 
empirical relations (cf. Appendix A) in order to compute the 
parameters of the pulsating star. The law of error propaga- 
tion will be employed to estimate the errors of the param- 
eters of the star: When one tries to estimate the value of a 
quantity z which is given as a function of other quantities 
z — f(xi, X2, ■ ■ ■) that can be measured directly, the error 
that should be attributed to the former one is given by 
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0.8 1.0 1.2 1.4 1.6 1.8 

M (in 1 .4 solar mass units) 



Figure 1. This diagram shows the range of mode-frequencies for 
twelve different equations of state, as functions of the normalized 
mass (M of the star. 



£ 



\dxi C 



dx i dxj 



-, 1/2 



(9) 



where a Xi is the measurement error of the quantity Xi and 
r Xi xj is the correlation between the errors of Xi and Xj. In 
Section 4 we will use the measured values of the frequen- 
cies of the / and the first p-mode (// and f p ) to obtain 
information about the parameters of the star. Since these 
are independent quantities (they are measured separately 
by different templates), the second sum of the above law 
will not be present in our calculations. The magnitude of 
the estimated errors of neutron star parameters will deter- 
mine the efficiency of the method and our chance to achieve 
our final goal, the determination of the equation of state of 
a neutron star. 



3 DETECTING PULSATION MODES 
3.1 Are the modes detectable? 

Two separate questions must be addressed in any discussion 
of gravitational- wave detection. The first one concerns iden- 
tifying a weak signal in a noisy detector, thus establishing 
the presence of a gravitational wave in the data. The second 
question regards extracting the detailed parameters of the 
signal, e.g., the frequency and e- folding time of a pulsation 
mode. To address either of these issues we need an estimate 
of the spectral noise density S n of the detector. 

The few pulsation modes of a neutron star that may 
be detectable through the associated gravitational waves all 
have rather high frequencies; typically of the order of several 
kHz. To illustrate this we show the mode-freque ncies f or all 
models considered by Andersson and Kokkotas (1998) as a 



function of the stellar mass in Figure [y. From this figure 
we immediately see that a detector must be sensitive to fre- 
quencies of the order of 8-12 kHz and above to observe most 
ui-modes. Of course, it is also clear that some equations of 
state yield ui-modes at somewhat lower frequencies. For ex- 
ample, for massive neutron stars with M ~ 1.8— 2. 3Mq (stiff 
EOS), as have been suggested for low-mass X-ray binaries, 
the u>-mode frequency could be as low as 6 kHz (see also re- 
cent results for the axial io-modes by Benhar et al (1999)). 



The p-modes lie mainly in the range 4-8 kHz, while all /- 
modes have frequencies lower than 4 kHz. This means that 
the mode-signals we consider lie in the regime where an in- 
terferometric detector is severely limited by the photon shot 
noise. For this reason a detection strategy based on resonant 
detectors (bars, spheres or eve n networks of small resonant 
detectors (Frasca & Papal995)) or laser interferometers op- 
crating in dual recycling mode seems the most promising. 
In fact, the range of mode-frequencies in Fig. ^ should moti- 
vate detailed studies of the prospects for construction dedi- 
cated ultrahigh frequency detectors. In the following we will 
compare three different detectors: The initial and advanced 
LIGO interferometers, for which 



J_ 



3/2 



Vf 



Hz 



-1/2 



(10) 



with h m = 3.1 x 1CT 22 , a = 1.4 and f m = 160 Hz for the 
initial configuration, and h m = 1.4 x 10~ 23 , a — 1.6 and 
fm = 68 Hz respectively , for the advanced configuration 
(Flanagan & Hughes 1995). We also consider an "ideal" de- 
tector that is tuned to the frequency of the mode and has 
sensitivity of the order of S n ~ 10 Hz ' (this is the 
sensitivity goal of detectors under construction ). It should 
be noted that the Advanced LIGO estimates are roughly 
valid also for spheric al de tectors such as TIG A, cf. Harry, 
Stevenson and Paik (1996). 



The detectability of the /, p and w-modes for differ- 
ent detectors can be assessed from (^|). The main problem 
in doing this is the lack of realistic simulations providing 
information about the level of excitation of various modes 
in an astrophysical situation. Still, given the frequency and 
damping rate of a specific mode we can ask what amount 
of energy must be channeled through the mode in order to 
be detectable by a given detector. We immediately find that 
detection of pulsating neutron stars from outside our own 
galaxy is very unlikely. Let us consider a "typical" stellar 
model for which the /-mode has parameters // = 2.2 kHz 
and Tf = 0.15 s. This corresponds to a IAMq neutron star 
according to the Bethe- John son eq uation of state, cf. Table 3 
of Andersson and Kokkotas ( 1995 ) . For this example we find 
that the /-mode in the Virgo cluster (at 15 Mpc) must carry 
an energy equivalent to more than 0.3Mqc 2 to lead to a 
signal-to-noise ratio of 10 in our ideal detector. Given that 
the total energy estimated to be radiated as gravitational 
waves in a supernova is at the level of 10 -5 — 10~ 6 Mqc 2 , we 
cannot realistically expect to observe mode-signals from far 
beyond our own galaxy. 

This means that the number of detectable events may 
be rather low. Certainly, one would not expect to see a su- 
pernova in our galaxy more often than once every thirty 
years, or so. Still, there are a large number of neutron stars 
in our galaxy, all of which may be be involved in dramatic 
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Table 1. The estimated energy (in units of A/qc 2 ) required in 
each mode in order to lead to a detection with signal-to-noise 
ratio of 10 from a pulsating neutron star at the center of our 
galaxy (10 kpc). 



Table 2. The relative errors in the extraction of the mode pa- 
rameters assuming a signal-to-noise ratio of 10. 



Detector 
Initial LIGO 
Adv LIGO 
Ideal 



/-mode 
4.9 X 10~ 5 
8.7 x 10~ 7 
1.4 x 10~ 7 



p-mode 
4.0 X 10~ 3 
7.0 x 10~ 5 
1.3 x 10~ 6 



w-mode 
6.8 X 10^2 
1.2 x 10" 3 
6.4 x 10" 6 



events (see the introduction for some possibilities) that lead 
to the excitation of pulsation modes. The energies required 
to make each mode detectable (with a signal-to-noise ratio 
of 10) from a source at the center of our galaxy (at 10 kpc) 
are listed in Table [l|. In the table we have used the data for 
the "typical" stellar model, for which the characteristics of 
the /-mode were given above, f p = 6 kHz and r v — 2 s, 
and f w — 11 kHz and t w = 0.02 ms. This data indicates 
that, even though the event that excites the modes must be 
violent, the energy required to make each mode detectable 
is not at all unrealistic. In fact, the energy levels required 
for both the /- and p-modes are such that detection of vio- 
lent events in the life of a neutron star should be possible, 
given the Advanced LIGO detectors (or alternatively spheres 
with the sensitivity proposed for TIGA). On the other hand, 
detection of ui-modes with the broad band configuration of 
LIGO seems unlikely. Detection of these modes, which would 
correspond to observing a uniquely relativistic phenomenon, 
requires dedicated high frequency detectors operating in the 
frequency range above 6 kHz. Still, we believe that the data 
in Table hi illustrates that neutron star pulsation modes may 
well be detectable from within our galaxy, and that the first 
detection may in fact come as soon as the first generation 
of LIGO detectors come on line. 





"■/// 


oy/t 


or (10 


/-mode 


8 x 10- 5 


0.2 


11.3 


p-mode 


3 x 10~ 6 


0.2 


3.8 


■iii-mode 


0.1 


0.3 


1.7 



IAMq stellar model we used in the previous section; see Ta- 
ble |^. We assume that each mode carries the energy required 
for it to be observed with signal-to-noise ratio of 10, cf. Ta- 
ble [!} (This is a convenient measure since it is independent 
of the particulars of the detector.) 

From the sample data in Table ^ one sees clearly that, 
while an extremely accurate determination of the frequen- 
cies of both the /- and the p-mode is possible, it would be 
much harder to infer their respective damping rates. It is 
also clear that an accurate determination of both the w- 
mode frequency and damping will be difficult. To illustrate 
this result in a different way, we can ask how much energy 
must be channeled through each mode in order to lead to 
a 1% relative error in the frequency or the damping rate, 
respectively. Let us call the corresponding energies Ef and 
E T . This measure will then be detector dependent, so we list 
the relevant estimates for the three detector configurations 
used in Table [j]. When the data is viewed in this way, cf. Ta- 
ble 3.2, we see that an accurate extraction of ui-mode data 
will not be possible unless a large amount of energy is re- 
leased through these modes. Furthermore, one would clearly 
need a detector that is sensitive at ultrahigh frequencies. In 
other words, it seems unlikely that we will be able to use the 
w-modes to infer the detail ed neutron star parameters as we 
have previously suggested (Andersson & Kokkotas 1996). 



3.2 How well can we determine the mode 
parameters? 

Let us now discuss the precision with which we can hope 
to infer the details of each pulsat ion m ode. After inserting 
Eq. (^) into formulae (B12) and (B16) we can compute the 
relative measurement error in the frequency and the damp- 
ing time of the waves by some appropriately designed detec- 
tor. After introducing a convenient parameter V according 
to 



V 



cl/2 
On 



-1/2 



10- 23 Hz-V2 J ^lOkpcy yiu -iviQ 
we find that the error estimates take the following form 



(11) 



£7 „ 
/ " 
and 



0.0042 V 



1 - 2Q 2 + < 



4Q 4 



lms/ 



— ~ 0.013 V 1 



10- 



Q- 



(lkHz) 



(12) 



(13) 



Also, for the time of arrival of the gravitational wave signal 
we get from ( B19 ) 
-l 



<jt ~ 0.0042 V~ ms 



(14) 



To illustrate these results we show the relative errors 
associated with the parameter extraction for the "typical" 



4 REVEALING THE EQUATION OF STATE 

In the previous section we discussed issues regarding the de- 
tectability of a mode-signal, and the accuracy with which 
the parameters of the mode could be inferred from noisy 
gravitational wave data. Let us now assume that we have 
detected the mode and extracted the relevant parameters. 
We then naturally want to constrain the supranuclear equa- 
tion of state by deducing the mass and the radius of the 
star (or combinations of them). In principle, the mass and 
the radius can be deduced from any two observables, cf., 
Table 2 of Andersson and Kokkotas (1998). In the absence 
of detector noise, several combinations look promising, but 
in reality only few combinations are likely to be useful. 

Consider the following example: We could in principle 
deduce the mass and the radius from a detected /-mode 
(assuming that we have extracted both its frequency and 
damping rate). However, as can be seen from Table 0, the 
estimated relative error in frequency is about three orders of 
magnitude smaller than the relative error in damping time. 
Hence, if these two measurements are to be used to deter- 
mine the mass and the radius of the pulsating star one must 
keep in mind that the measurement of frequency is far more 
accurate than the measurement of damping time. And it is 
clear that by combining // with r/ we will only get accurate 
estimates for M and R if the energy in the mode is sub- 
stantial. The same is true for the combination f p and t p , as 
well as any combination involving the w-mode data. Still, we 
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Table 3. The estimated energy (in units of Mqc 2 ) required in each mode in order to lead to a relative error 
of 1% in the inferred mode-frequency (Ef) and damping rate (E T ). In cases where no entry is given, the 
required energy is unrealistically high (typically larger than Mqc 2 ). The distance to the source is assumed 
to be 10 kpc. 

Detector /-mode p-mode w-modc 





Ef 




E T 




Ef 




E T 


Ef 


E T 


LIGO 1 


3.1 x 10" 


9 


1.9 x 10" 


-2 


2.8 x 10" 


10 








Adv LIGO 


5.5 x 10" 


11 


3.5 x 10- 
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5.6 x 10" 
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8.8 x 10" 
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9 x 10~ 4 


6.1 x 10~ 3 



should not discard the possibility that there will be unique 
events for which w-modes ca rry the bulk of the ene rgy, as 
for example the scattering of (rbminariga et al. 199S) or the 
infall of a smaller mass on neutron stars ( porclli 1997| ). In 
such cases , the strategy proposed by Andersson and Kokko- 
tas (|l998|) will be useful. 

From the data in Tables ^j.2[ it seems natural to use 
the p-mode in any scheme for deducing the stellar param- 
eters. But, as we have already mentioned, there does not 
exist a "nice" relation between frequency and stellar param- 
eters for the p-mode, that is independent of the equation of 
state. Still, the strategy that would seem the most promis- 
ing on the basis of the present results is based on using 
the frequencies of /- and p-modes. A possible method is 
based on the following steps: The first step is to invert the 
empirical relation for the /-mode frequency, Eq. (Al), in 
order to transform the measurement of ff to an estimate 
of the mean density of the star. Let us define a parameter 
x = (M/R 3 ) 1 ' 2 , where M = M/1AM® and R = R/10 km, 
in order to simplify the notation. Then, once // has been 
measured x can be computed from 



_ //(kHz) -0.78 
X ~ L63 ' 
with a corresponding relative error 

0.97 1ms 1 



(15) 



r f ff- 0.78 



+ 



0.01 



ff - 0.78 



+ 0.006^ 



1/2 



where ff is expressed in kHz. 

The first of the three terms in the square root comes, ob- 
viously, from the measurement error of ff (the extra compli- 
cating factor, related with the Q of the /-mode, that arises 
there when one tries to express the result with respect to 
the signal-to-noise, has been simplified to unity since the Q 
of the /-mode is generally a large number). The other two 
terms arise from dispersion of the data for the various equa- 
tion of state. For a typical /-mode frequency of 2 kHz, and 
damping time 0.2 s the relative error in x is ~ 0.01 (assum- 
ing p > 6). Actually, the first term in Eq. ( |l6[ ) is negligible 
and could be omitted. Therefore 



(17) 




+ 0.006 2 



Here, we need not worry about the sign of the factor 
ff — 0.78; it is always positive since // > 1.4 kHz for every 
ste llar model in the dataset, see Figure [I] here, or Figure 1 
of (Andersson & Kokkotas 199£ ) . 



Then, by measuring the frequency of the first p-mode 
— which is expected to c arry roughly as much energy as the 
/-mode, see Allen et al (1998), one could place an error box 



i — 1 — i — 1 — r- 



~ i — 1 — i — 1 — i — 1 — i — 1 — i — 1 — i — ■ — r 




Error Box 



0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 
(M/R 3 )" 2 



Figure 2. This diagram shows how the observable frequency of 
the first p-mode scales with the square root of the mean density of 
a neutron star, (M/R 3 ), according to numerical data for various 
equations of state. 



(16) 

on a diagram of f v vs. x, where all theoretical models for 
the equation of state are drawn. We illustrate this in Fig- 
ure ^. This way we can identify the most likely equation of 
state. Detecting gravitational waves from pulsating neutron 
stars ensures quite accurate measurements of // and f p , so 
the error box in Figure |^ is remarkably small. Hence, this 
method can easily distinguish between the different equa- 
tions of state in our dataset. In fact, what makes this method 
so efficient is the fact that different equations of state are 
described by quite distinct curves in an f p — x diagram. 

Finally, we would like to infer the mass and radius of 
the neutron star. To do this we can use the data in Fig- 
ure |^, which gives the connection between the p-mode fre- 
quencies of the different equations of state and the stellar 
compactness. From this diagram, we can immediately use 
the detected f p to infer the compactness of the star, once 
the right equation of state has been identified. Having esti- 
mated both the average density and the compactness, it is 
an elementary calculation to obtain the mass and the radius 
of the star. 

At this point it is natural to ask the following question: 
What if the true equation of state is not close to one of the 
present models? This may well be the case. After all, our 
understanding of the state of nuclear matter at supranu- 
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Figure 3. This diagram shows how the observable frequency 
of the first p-mode scales with the compactness of a neutron 
star, (M/R), according to numerical data for various equations 
of state. 



clear densities still awaits observational testing. Should the 
equation of state be markedly different from the ones in our 
sample we will get an error box which does not lie close to 
any of the equations of state in Figure ^| Should one then 
want to estimate the stellar parameters, one could construct 
a sequence of polytropes (p = Kp 1 ,7=1 + 1 /n) for various 



values of K and n (see for example the graphs in (Ander 
sson fc Kokkotas 199?j )). The appropriate combination of 
these free parameters will bring the corresponding equation 
of state within our error box. Then, one could use Figure 
and read off the "correct" compactness of the star (M/R) 
and from this compute its radius and mass. 



5 DETERMINING THE POSITION OF THE 
SOURCE 

As with other kinds of gravitational- wave sources, a network 
of at least three detectors is needed to pinpoint the location 
of the source in the sky. The difference in arrival time for 
the three detectors could be used to determine the position 
of the source. The higher the accuracy in measuring the 
time of arrival at each detector, the more precise will be 
the positioning of the source. Two remote detectors, at a 
distance d apart from each other will receive the signal with 
a temporal difference of 



AT = - cos ( 

c 



(18) 



where c is the speed of light, and 9 is the angle between the 
line joining the two detectors and the line of sight of the 
source. Therefore, the accuracy by which this angle can be 
measured is 



The y/2 arises from the measurement errors of the two times 
of arrival. If one assumes an 'L' shaped network of 3 detec- 
tors with arm length of d = 10, 000 km, Eqs. @ and @ 
lead to an error box on the sky with angular sides of 1°, 
at most (for specific areas of the sky, and large signal-to- 
noise ratios the angular sizes could be much smaller). This 
is quite interesting since one could then correlate the detec- 
tion of gravitational waves with radio- , X-ray or gamma-ray 
observations directed on that specific corner of the sky. 



AO = 



\/2' 



<ttc 



dsin 8 



(19) 



CONCLUDING REMARKS 

In this paper we have extended previous studies of the de- 
tectability of gravitational waves from pulsating neutron 
stars to include the statistical errors associated with an 
analysis of a weak signal in a noisy data stream. We have 
shown that the generation of detectors that is presently un- 
der construction may well be able to observe such sources 
from within our own galaxy. Detections from distant galax- 
ies seem unlikely, unless a sizeable amount of energy (of the 
order of Q.QIMqc?) is released through the pulsation modes. 
This means that we expect the event rate to be rather low. 
One would certainly not expect to see more than perhaps 
three neutron stars being born in supernovae per century 
in our galaxy. Of course, other unique events in a neutron 
stars life, like starquakes and phase-transitions, may lead to 
relevant signals. Perhaps the most interesting possibility is 
the possible association between gravitational waves from 
pulsation modes and gamma rays from the soft-gamma re- 
peaters (the magnetars). The number of such events is very 
hard to estimate at the present time. Given the many un- 
certainties in the available models of all these scenarios, we 
certainly cannot rule out the possibility of detecting gravi- 
tational waves from pulsations in neutron stars. 

However, and this is an important point, the chances of 
detecting pulsating neutron stars would be much enhanced if 
one could construct dedicated high-frequency gravitational- 
wave detectors sensitive in the range 5-10 kHz. This provides 
a serious challenge to experimenters, but that the pay-off 
of a successful detection of oscillating neutron stars would 
be great is clear from our present analysis of the inverse 
problem. We have shown that the oscillation frequencies of 
the fluid / and p modes can be accurately deduced from 
detected signals. Other parameters, like the particulars of 
the gravitational-wave w modes will be much more difficult 
to infer, unless the signal-to-noise ratio of the detection is 
unexpectedly large. This means that a previously proposed 
strategy for deducing the parameters of the star (the mass 
and the radius) from observations of the / and a ui-mode is 
unlikely to be practical. However, we show that an equally 
powerful strategy can be based on detected / and p-modes. 
Given an observation of / and p-mode oscillation frequen- 
cies, with the estimated accuracies, we can easily rule out 
many proposed equations of state. In other words, we have 
proposed a gravitational- wave "fingerprint analysis" for neu- 
tron stars holds a lot of promise, and may help us reveal the 
true equation of state at supranuclear densities once gravi- 
tational waves are observed. 
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APPENDIX A: EMPIRICAL RELATIONS 

In this Appendix we list the empirical relations deduced 
from data for twelve different realistic equations of star. 
These relations between the observables (frequency and 
damping of the modes) and the stellar parameters (mass 
M and radius R) are essenti ally th e same as the ones listed 
by Andersson and Kokkotas (1998), but here we also include 
the relevant statistical errors. 

We have constructed empirical relations for both the 
fundamental mode of fluid pulsation (the /-mode) and the 
slowest damped of the gravitational wave ui-modes. For the 
/-mode, the frequency scales with the mean density of the 
neutron star according to 



(0.78 ±0.01) + (1.63 ± 0.01) (i| 



1/2 



(Al) 



while the damping time of the /-mode can be described by 



Is 



M 3 
W 



(22.85 ± 1.51) - (14.65 ± 1.32) { M- 

R 



,(A2) 



where M = M/1AM Q and R = R/W km. The indicated 
uncertainties are simply the statistical errors that show up 
if one takes into account all stellar models in th e da ta set, 
on equal basis. In principle, equations (Al) and (A2) could 



be inverted to compute the average density (M/R 3 ) and 
compactness (M/R) of the star from measured values of 
// a nd T f. However, this procedure proves unreliable since 
Eq. (A2) is a double- valued function with respect to com- 



pactness, and the presence of errors makes it impossible to 
infer the compactness of the star with reasonable accuracy, 
even if the characteristics of the waves could be measured 
with extreme precision. 

The w-modes are pulsations directly associated with 
spacetime itself. They have relatively high oscillation fre- 
quencies (6-14 kHz for typical neutron stars) and barely ex- 
cite any fluid motion. They are also rapidly damped, with 
typical lifetime of fraction of a millisecond. The frequency 
of the first w-mode scales with the stellar compactness as 



1kHz R 



(20.95 ±0.33) - (9.17 ±0.29) 



(A3) 



while the damping rate of the mode is well described by 



M 



, M 

(3.90 ± 4.39) + (104.06 ± 8.33) 

Ft 



(67.28 ± 3.84) 



(A4) 



APPENDIX B: THE FISHER AND 
COVARIANCE MATRIX 

For a typical mode signal the components of the symmetric 
Fisher matrix, defined by 

dh , dh 



where 9% = (e, rj, £, £) are 



r ££ = 



r - 



1 + 24Q 4 + 32Q 6 2 
(1±4Q 2 ) 2 9 

3-4Q 2 2 
2(1 + 4Q2)2 P ' 



T £C = -2QV , 
r 1 2 



1+4Q 2 
3 + 6Q 2 + 8Q 4 2 



(1 + 4Q 2 ) 2 



r l r? 



3 + 4Q 2 2 
2(1 + 4Q 2 ) P 



(Bl) 

(B2) 

(B3) 
(B4) 
(B5) 

(B6) 
(B7) 
(B8) 
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r cc = (l+4QV , (B9) 

r« = o , (bio) 
r« = p 2 ■ (Bii) 

The inverse of the Fisher matrix, E;j = r~- , the so- 
called covariance matrix, has the following components: 

_ 1-2Q 2 +8Q 4 1 
2Q4(1 + 4Q2) P 2 ■ 

3 - 4Q 2 1 

Q2(1+4Q 2 )^2 

-1+4Q 2 1 
jeC ~ 2Q2(l + 4g2)p2 

-1 + Q 2 1 
2Q 4 p 2 



(B14) 



E«e = -^TJ . (B15) 



_ 4(5 + 4Q 2 ) . 



^=(TW)7' (B17) 

-3-2Q 2 1 



(1 + 4Q2 


) p 2 


-4 


1 


(1+4Q2) 


P 2 : 


-3 - 2Q 2 


1 


Q 2 


P 2 ' 


2 


1 


(1 + 4Q2) 


P 2 ' 


1 1 








(1 + 2Q 2 ) 2 1 


2Q 4 


P 2 



N - ; ,, - ( B1 9) 

(B20) 

1 ' I •'" ! : moil 

« - 2Q 4 ^ ' (B21) 

These relations can be used to estimate the measurement 
error in the various parameters of the mode-signal; see dis- 
cussion in the main text. 
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